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Low-density nuclear matter at finite temperature is considered representing the strong coupling 
situation of a highly correlated fermion system. One-particle self-energies and the density of states 
in the vicinity of the pairing transition point are presented. Within the Green function approach 
model calculations are performed using the thermodynamic T matrix in ladder approximation. As 
a main result the formation of a pseudogap in the level density near the critical temperature of the 
pairing phase transition has been found. It is shown that compared to mean field, the picture of the 
onset of the pairing transition is essentially changed if correlations are taken into account. 
Keywords: pseudogap, level density, pairing, superfluidity, strongly coupled fermion systems, nuclear 
matter. 
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Among strongly coupled fermion systems which show 
a transition to a superfluid or superconducting state 
|J such as, e.g., the electron gas in metals, condensed 
He , or electron-hole pairs in excited semi-conductors, 
low-density nuclear or neutron matter turns out to be 
a particularly strongly correlated system. According 
to BCS theory {§ the ratio A/e F (A being the BCS 
gap, ep the Fermi energy) is of the order 1CP 3 in or- 
dinary metals whereas it is 0.1 in high T c superconduc- 
tors. In finite nuclei the ratio is about 0.03 and thus 
much closer to the scenario of high T c superconductors 
than to normal superconductivity. Indeed, the coher- 
ence length £ = Upft/A for nucleons in the spin sin- 
glet (S — 0), isospin triplet (T — 1) channel becomes 
of the same order as the internucleon distance at densi- 
ties p « po/10. . . po/5 with po = 0.17 fm _ the nuclear 
matter saturation density. In the S = 1, T = neutron- 
proton channel (the deuteron), where the interaction is 
more attractive, the correlations are even stronger. 

Referring to this, a very fundamental quantity is the 
two-particle spectral function which is governed by such 
correlations. At the same time, the two-particle spec- 
tral function is an important ingredient to determine 
single-particle properties. In particular, in high-T c mate- 
rials such as underdoped cuprates, strong correlations are 
present in the electron-hole system of the normal (non- 
superconducting) phase leading to the formation of the 
so-called pseudogap in the charge and spin excitation 
spectra. It is known from experiments and from theo- 
retical studies that the pseudogap occurs below a certain 
temperature T* having similar characteristics as the su- 
perconducting gap which grows continuously out of the 
pseudogap ||. It appears that the pseudogap is a pre- 
cursor of the superconducting gap. 

In this letter we focus on low-density nuclear matter 
at finite temperatures. Our main point of interest is 



the isospin singlet pairing |4],|5| being of relevance mostly 
in symmetric nuclear matter 1| which will be studied 
here. It should be pointed out that low-density nuclear 
or neutron matter not only exists during stellar collapse 
or in the outer sphere of neutron stars but is an impor- 
tant prerequisite to study pairing effects via the rather 
widespread and well accepted local density approxima- 
tion also in the surface layer of real nuclei. There, evi- 
dence of strong correlations is given by electron scatter- 
ing experiments off nuclei M where the nucleon spectral 
function and momentum distributions are determined. In 
such highly correlated systems the mean-field BCS theory 
is not the appropriate tool for the description of pairing 
effects. 

It is therefore of great importance to include correla- 
tions in the theoretical treatment of such systems. We 
use the T matrix approach within the Green function the- 
ory (III . This is a very powerful tool to describe bound 
states (deuterons) in the normal state and the pairing 
transition (formation of Cooper pairs) ||,[Hj-[l3| . 

The single-particle propagator obeys the Dyson equa- 
tion fii 



Gi(0,z) = 



Opp> 



z — p 2 /(2m) — S(p, z) 



(i) 



where the self-energy £ contains the influence of the 
medium on the propagation of a single particle which is 
systematically accounted for in the single-particle spec- 
tral function A(j>,uj), based on which the level density is 
given by 



W(w) = £>(p,< 



(2) 



In the following we are interested in energies around the 
Fermi level and therefore we will neglect the imaginary 



part of E, i.e. the spectral function is approximated by 
a (5-function at the quasi-particle energy e(p) given by 



e(p) = ^ + i7(p) ; J7(p) = Re E(p, e(p)) . (3) 
Then the level density can be calculated via 

JV(w)=5>(w-e(p)), ( 4 ) 

p 

which is identical to (||), see e.g. [1J]. The calculation of 
(0) yields the final expression for the density of states 



N(u) 



1 de (p) 
p 2 dp 



e(p)- 



(5) 



It is completely determined by the real part of the on- 
shell self-energy and related to the effective mass m* (p) = 
2 N{ijj)/p. The consistent inclusion of Im £ is in progress. 
Within a cluster decomposition |l5| the single-particle 
self-energy can be calculated from the two-particle T ma- 
trix 72(12, 1'2', z) taking only two-particle correlations 
into account. One obtains for the real part pi 

Re£(pi,w) = Rc£ (1) (p 1 ,w)+ReI] (2) (> 1 ,Lj), (6a) 

ReS( 1 )(p 1 ,^) = ^/(2)ReT 2 (12,12,^ + 6(p 2 )), (6b) 

Re #*> fo, a,) Jyf^L ^)^ a (12, 12, a/ + iO) 

with /(2) = /(e(p 2 )) - [exp{(e( P2 ) - M )/T} + l]" 1 being 
the Fermi and #(£) = [exp{(£ - 2\i)jT} - l]" 1 the Bose 
distribution function. The two-particle ladder T matrix 
is given by the Galitskii-Feynman equation B 

T 2 (12, 34, z) = V(12, 34) + ^ V(12, 1'2') 



1'2' 



xG°(l'2',z)r 2 (l'2',34,z). 



(7) 



The uncorrelated two-particle Green function G 2 is given 
as a product of two single-particle propagators via the 
one-particle spectral function. Using the quasi-particle 
approximation as mentioned above one obtains 



G° 2 (12,z) 



1 - /(l) - /(2) 
z - e(pi) - e(p 2 ) 



(8) 



Equations (0) and (||) together with (0) and rf|) form a 
self-consistent set of equations which is solved iteratively. 
The usual first-order Briickncr-Hartree-Fock approxima- 
tion M] can be recovered replacing in (|8|) the expression 
1 - /(l) - /(2) by (1 - /(1))(1 - /(2)), i.e., neglecting 
the "backward going" part. As we will demonstrate later 
the inclusion of the hole-hole sector is, however, essential 
to account for precritical pair fluctuations. 

For an exploratory calculation we use a separable rep- 
resentation of the nucleon-nucleon potential according to 
Yamaguchi fll7| 



V a (p,p) = —£v(p)v(p') ; v(p) 



1 



p 2 + 13 2 



(9) 



with the relative momenta p and p' , a denotes the inter- 
action channel, A Q the coupling strength of the respective 
channel, and Q the normalization volume. The param- 
eters = 1.4488 fin -1 and \ a = 151.7 MeVfm 3 for the 
1 So channel and A Q = 216.1 MeVfm 3 for the 3 Si channel 
are chosen to reproduce the low-energy scattering phase 
shifts as well as the deuteron binding energy. This po- 
tential yields somewhat too strong pairing properties as 
compared to those calculated, e.g. with the Paris poten- 
tial |5J , but this will not invalidate the important conclu- 
sions which we will draw from our study below. Using 
the potential (||) a solution of (0) is particularly simple. 
Transforming to relative p = (p{ — p 2 )/2 and total mo- 
mentum P = p{ + P2 one obtains 



T a (pp'P,z) 



V a (p,p') 



(10) 



1 - J a (P,(i,T,z) ' 
J a {^T,P,z) = ^ fdpp 2 v 2 (p)G a 2 ( P ,P,z). (11) 



Since Re £( 2 ) ( |6q) is of particular interest we will con- 
sider it in more detail. In particular, there is a connection 
to the two-particle spectral function A 2 that is defined 
as the imaginary part of the full two-particle Green func- 
tion G 2 . Instead of A 2 {pip 2 ,p' 1 p' 2l Lo) we will consider the 
contraction (pip = P/2 ±p, p[ 2 = P/2 ±p') 

F 2 (P,u) = ^ijp)l^G 2 {p 1 p 2 p'iP' 2 ^ + ^)v{p'), (12) 

which is a function of only the total momentum and the 
energy. With (12) one can rewrite (pq) 



Re£( 2 >( Pl , W ) = W ^U-i(l2,12)4 



g (q/ )F 2 (P,q/) 
lu - e{p 2 ) 

(13) 



It can be shown that F 2 obeys the following energy- 
weighted sum rule (the analogue to the /-sum rule in 
the particle-hole channel) 



+00 



p 2 



m 4m 



Q(p,P)v 2 (p) 



+ 



d 3 P 2 

(2tt) 3 



v 2 (p)Q(p,P)) (14) 



(Q denotes the angle averaged Pauli blocking operator). 
In the following we are going to present results of a 
calculation for symmetric nuclear matter at the density 
p = po/3. Using the self-consistent procedure described 
above one obtains from the Thouless criterion of the T 2 
matrix pd a critical temperature T c = 3.6 MeV in the 
3 Si-channel. 
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FIG. 1. Single-particle potential (g) of the nucleon in nu- 
clear matter as a function of the momentum p over the Fermi 
momentum pf — 0.92 fm _1 at the density p = po/3 and tem- 
perature T = 5 MeV. The contributions LA 1 ' and U ^ are 
defined according to Eqs. (M), (pbh, and (pq), respectively. 

In our example we chose the temperature T = 5 MeV 
being close but above this critical value. In Fig. [I] the 
self-consistent real part of the on-shell self-energy (single- 
particle potential) U{p) is displayed as a function of the 
momentum p showing the contributions U^ and U^ 
and their total sum according to Eqs. (pq), (pq), and 
(paj), respectively. The sum U(p) — U^(p) + U^(p) is a 
delicate balance of the first and the second contribution 
to the single-particle potential. Although each contribu- 
tion shows a smooth behavior the sum of both exhibits 
two plateaus around the Fermi momentum pp sandwich- 
ing a steeper slope. This effect can be traced back to 
the pairing instability of the T 2 matrix that is very pro- 
nounced already close to T c (see pT]Jl3| ) ■ A similar result 
has been found at zero temperature by Vonderfecht et al. 
|L8| who calculated U(p) in the same approach but un- 
der consideration of the BCS gap in order to avoid the 
pairing instability. 

Because of the subtlety of the effect it is important to 
check numerical accuracy. We therefore calculated the 
real part of the self-energy also via dispersion relation 
from the imaginary part according to 



ReS(p 1 ,e(p 1 )) = S HF (p 1 )+|^ 



cIlu Im £(pi , u> + iO) 



(15) 



u - e(pi) 
S HF denotes the Hartree-Fock shift and Im £ is given by 

/j3 
-^ImT 2 (12, 12, lo + e{p 2 ) + »0) 

x[/(e(pa))+s(« + c(p2))]. (16) 

Within the bounds of numerical accuracy the calculation 
of Re £ according to (|l5|) , (jl^) yielded equal results com- 
pared to the procedure following from (£ 
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FIG. 2. Nucleon level density (J5j) as a function of the en- 
ergy u> — H normalized to the chemical potential fi given for 
the density p = po/3 and temperature T = 5 MeV. Re- 
sults of three different approximations are shown: T matrix 
approach (solid) and for comparison Briickner-Hartree-Fock 
(BHF, dashed) and mean-field BCS theory (dotted). The 
dot-dashed curve is a T = calculation in BCS approxima- 
tion. 

According to (|5|) the behavior of the single-particle po- 
tential is closely related to the density of states. In Fig. g 
N(lu) is plotted as a function of u> — fi for the same val- 
ues of density and temperature as in Fig. |l|. One observes 
a depletion of the density of states around the chemical 
potential and the rise of two pronounced peaks. This for- 
mation of a so-called pseudogap (see below) in the level 
density is a precursor signal of the superfluid phase tran- 
sition expected to occur at T c — 3.6 MeV. For compari- 
son, in Fig. results of other approximations are given. 
The mean-field BCS approach at T = 5 MeV (dotted 
curve) which is closely below the critical temperature of 
T C BCS = 5.08 MeV exhibits a tiny gap of twice the BCS 
gap at Fermi momentum [A(pp) = 2.1 MeV] around the 
chemical potential. The BCS result for zero temperature 
(dot-dashed) shows a considerably larger gap. There is a 
striking similarity between the BCS result at zero tem- 
perature and the T matrix pseudogap with respect to the 
location of the two peaks. We also calculated the level 
density in Briickner-Hartree-Fock (BHF) approximation 
|16[ which neglects the hole-hole contribution to the ef- 
fective interaction which is necessary for the description 
of pairing. Interestingly, the BHF calculation yields an 
opposite behavior compared to the T matrix approach 
which is known as the enhancement of the effective mass 
at the Fermi surface pqj . 

In the following we want to discuss the consequences 
of our calculation in detail. The occurrence of the pseu- 
dogap in the density of state is a phenomenon which has 
recently been discovered experimentally and discussed 
theoretically in the context of high T c superconductors 
B. Some authors attribute the measured depression in 



the density of states to pair fluctuations in analogy to 
the way described here Q. In the nuclear physics con- 
text our finding is novel and certainly opens up a so far 
unrevealed aspect of the nuclear density of states. 

However, the situation in finite nuclei is more com- 
plex. One should realize that all the effort over the last 
decades concerning the nuclear level density parameter 
rather went into the opposite direction trying to explain 
why the experimental single-particle level density of finite 
normal fluid nuclei around the Fermi energy is enhanced 
rather than diminished with respect to the mean-field 
result. Indeed, in including long range particle-hole cor- 
relations such as surface vibrations, it is very well known 
p0[ that the single-particle potential is flattened in the 
vicinity of the Fermi energy rather than steepened as 
here. The strong collectivity of surface vibrations in nu- 
clei is certainly a quite particular feature of the finiteness 
of the system which makes that finally the particle-hole 
correlations in nuclei win over the particle-particle ones 
(otherwise there would be no deformed nuclei). However, 
the location of the transition point to stable deformation 
results from a detailed balance between the battle of pair 
and shape fluctuations. We believe that apart from rare 
occasions [pll the study of the influence of pair vibra- 
tions on the nuclear level density is very much under- 
represented in the nuclear physics literature and further 
investigations are of great importance for the detailed 
comprehension of nuclear systems. Indeed, e.g. for finite 
nuclei at T = the average level spacing takes over the 
role of temperature. There exists a critical level spacing 
from where on the nucleus becomes superfluid. Before, 
there exist precritical pair fluctuations quite analogous 
to the situation in infinite matter at T ~ T c . On the 
other hand in extended objects like neutron star matter 
it is not clear whether particle-hole correlations win over 
particle-particle correlations at all. The widely believed 
scenario where the effective nucleon mass in nuclear mat- 
ter is enhanced at the Fermi surface |lq | needs a more 
critical analysis once one allows for pairing correlations. 
This is a very important issue for future studies. 

In short, in this work we have for the first time per- 
formed a self-consistent Galitskii-Feynman-Hartree-Fock 
calculation at finite temperature approaching the super- 
fluid phase transition point T c from above. In parallel 
with theoretical work concerning high T c superconduc- 
tivity we also found in low-density nuclear matter the 
formation of a pseudogap in the density of states. This 
new finding may have not only important consequences 
for the physics of neutron and collapsing stars but may 
also effect details of finite nuclei. Further work is needed 
to settle the subtle competition between pair and den- 
sity fluctuations in this respect. Our results open up an 
essential new aspect in the important debate around the 
nuclear density of states. 

We acknowledge useful discussions with M. Baldo, W. 
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Preliminary results with respect to the pseudogap in nu- 
clear matter have been shown in [E2fl. 
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